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Abstract 

A BPS exact wall solution is found in M = 1 supergravity in four 
dimensions. The model uses chiral scalar field with a periodic superpo- 
tential admitting winding numbers. Maintaining the periodicity in super- 
gravity requires a gravitational correction to superpotential which allows 
the exact solution. By introducing boundary cosmological constants, we 
construct non-BPS multi-wall solutions for which a systematic analytic 
approximation is worked out for small gravitational coupling. 
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Introduction 

Supersymmetry (SUSY) has been most promising for unified theories beyond the stan- 
dard model |]. Recently new possibilities have been added by the "Brane World" sce- 
nario, where our four-dimensional spacetime is realized on topological defects such as 
walls embedded in a higher dimensional spacetime 0. More intriguing model has also 
been proposed by using the so-called warped metric in a bulk AdS spacetime of higher 
dimensional gravity together with an orbifold S 1 /Z 2 with a boundary cosmological con- 
stant finely tuned with the bulk cosmological constant ||. The warped spacetime has 
provided an attractive notion of localized graviton on the wall [|]]. Much efforts have been 
devoted to implement supergravity in orbifolds |J, work out solutions in gravity with 
a scalar field to replace the orbifold by a smooth wall configuration in four dimensions 
H [7j and in higher dimensions ||7|-^2|. It has been noted that the set of second order 
differential equations for coupled scalar and gravity theories can be reduced to a set of 
first order differential equations using a potential similar to the superpotential even for 
nonsupersymmetric case ||. 

In SUSY theories, single wall can preserve half of the SUSY, and is called a BPS state 
Hy| . Properties of topological defects such as walls have been extensively studied in SUSY 
theories fI3fl- |18|| and exact solutions have been useful as illustrated by a 1/4-BPS junction 



of walls lis] . Coexistence of BPS walls conserving different supercharges can break SUSY 
completely . When SUSY is broken, the stability of the multi-wall system is no longer 



guaranteed. We have found that topological quantum number such as winding number can 
stabilize the non-BPS multi-wall configurations. In an M = 1 globally SUSY sine-Gordon 
model, we have constructed an exact non-BPS multi-wall solution with a nonvanishing 
winding number which has no tachyon and is stable in a compact base space with the 
radius R |T7]]. We have also found an example of a stable non-BPS bound state of BPS 



walls in an extended model Since the compact space radius R becomes a dynamical 
quantity in gravity theories, we are led to ask if the stabilization mechanism due to the 
winding number still works in supergravity. To study properties of walls in supergravity 
theories, we first need to find solutions of BPS and non-BPS walls. 

The purpose of our paper is to present an exact BPS solution of the SUSY sine-Gordon 
model embedded into M = 1 supergravity in four dimensions and a systematic analytic 
approximation for non-BPS multi-wall configurations in powers of gravitational coupling. 
We find that a gravitational correction to the superpotential is needed to maintain the 
periodicity of our SUSY sine-Gordon model as an isometry realized by a Kahler transfor- 
mation. This allows us to obtain an exact BPS wall solution in the Af = 1 supergravity. 
Since equations of motion for matter and gravity are simultaneously satisfied, our solution 
requires no artificial fine-tuning of parameters such as boundary and bulk cosmological 
constants contrary to the original orbifold model 0. We discuss classification and global 
properties of non-BPS solutions using the first order differential equations of Ref.||, and 
succeed in constructing multi-wall configurations with and without winding number us- 
ing scalar field. We show that multi-wall configurations require negative energy density. 
We introduce a negative energy density as a boundary cosmological constant similarly 
to p, U). The scalar field satisfies the equation of motion without any external sources 
and is smooth everywhere. We also obtain a systematic analytic approximation in powers 
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of gravitational coupling and express necessary boundary cosmological constant in terms 
of elliptic integrals, for instance. To this end, we identify a set of first order differential 
equations for walls with single scalar field corresponding to the no gravity limit of the 
method of ||. We analyze the region of large gravitational coupling numerically and find 
qualitatively similar results. 

Exact BPS Wall Solutions with a Periodic Potential 

To avoid inessential complications, we will consider three-dimensional walls in four- 
dimensional theory to obtain BPS and non-BPS solutions except stated otherwise. Pre- 
viously we used a simple model admitting the winding number in M = 1 SUSY theory, 
the SUSY sine-Gordon model [[H]] with the Kahler potential K for the minimal kinetic 
term and the periodic superpotential P for a chiral scalar superfield $ = ((f), if), F) 



K(&, $) = $t $j P global ($) = ^ sin (|$) . 
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Coupling the model with the minimal kinetic term to M = 1 supergravity gives a 
Lagrangian whose bosonic part is given in terms of the superpotential P(4>) with scalar 
field replacing the superfield $ 
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where g mn is a metric of the spacetime, g = detg mn , R is the scalar curvature, and M P is 
the Planck massQ . 

Since we wish to study the winding number of 0, the periodicity under (|2|) should be 
maintained. The Kahler potential transforms under (EJ) as 



fl- \ 9 J 



(4) 



The Lagrangian is invariant under such Kahler transformations in the globally SUSY the- 
ory. In supergravity, however, it has been known that the superpotential should transform 



in a definite way under Kahler transformations to make the theory invariant |J|] . We find 
the necessary gravitational corrections to the superpotential which realizes the desired 
periodicity (fj) as 



K(&, $) = P($) = e ^P g lobal($) 
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x We follow the convention of Ref.[[l9| except that the Einstein indices are denoted by m, n, . 
dimensions, and by /j,, v, . . . in three dimensions. 



(5) 

in four 
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This gravitational correction is crucial to obtain exact periodicity and also an exact solu- 
tion of BPS walls. 

We find that the SUSY vacuum condition in supergravity gives precisely the same 
SUSY vacua as in the globally SUSY case, = g</>/A = f + nir, nfZ 
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No other vacua are allowed in spite of nontrivial dependence on Mp. By requiring half 
of the SUSY to be conserved, one obtains the BPS equations in Af = 1 supergravity We 
assume the following warped metric Ansatz || 



ds 2 = g mn dx m dx n = e 2A{y) r]^dx^dx v + dy 2 , (//, u = 0, 1, 3) 



(7) 



where A(y) is the warp factor, r]^ u = diag.(— , +, +), and the extra dimension is denoted 
as x 2 = y. Parametrizing the conserved SUSY parameter as ((y) = e ld ^a 2 ^(y), the half 
SUSY condition for gravitino gives the BPS equation for the Killing spinor ( a and the 
warp factor A [BJ 
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The BPS equation for the scalar field reads 
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Let us solve these BPS equations from y = — oo choosing a SUSY vacuum <fi = —n/2 
and 9 = n/2 as the initial condition at y = — oo. We observe that the right-hand side of 
the BPS equation for the imaginary part ((f) — 4>*)/2i is proportional to the imaginary part 
itself if 9 = it/2. Then the right-hand side of the BPS equation fliCf ) for the imaginary part 
(0 — 4>*)/2i vanishes at y — — oo. Similarly the right-hand side of (0) for the phase 9(y) 
vanishes at y — — oo. Therefore these first order differential equations dictate vanishing 
imaginary part of (ft and the constant phase of ( a for any y : 9(y) = ir/2, <p(y) = 4>(y)*. 
Then equation (|8|) can be integrated to give the Killing spinor ( a = — = ie^ A . The 
remaining BPS equation for = 0* becomes identical to the globally SUSY case and that 
for A becomes 



d<f) A 2 g dA 
dy g A dy 
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Anti-BPS equation conserving the opposite SUSY charges is given by 9 = —tt/2. We 
find solutions for a wall localized at y$ connecting from = —tt/2 + nir at y = — oo to 
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= 7r/2 + mraty = oo for any n € Z as (anti-)BPS solution for n = even (odd) as shown 
in Fig.[I] 
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where we suppressed an integration constant for A which amounts to an irrelevant nor- 
malization constant of metric. The Killing spinor is given with the normalization N 
( a = iN [cosh A(y — yo)]~ A ^ 2g Mp ^ '■ We have verified that these BPS solutions indeed 
satisfy the equations of motion. The energy density of the scalar field of the (anti-)BPS 
wall is given by 4T 3 = 4A 3 /g 2 . 
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Figure 1: BPS wall solution <p(y) and A{y) connecting 



-tt/2 to 4> = tt/2. 



Let us examine the thin- wall limit ]7||. By fixing the energy density 4T 3 and taking 
a limit g — ► oo, we recover the Randall- Sundrum model |4j]. The kinetic energy of the 
scalar field \d(f>/dy\ 2 reduces to the half of the boundary cosmological constant Ab OU ndar y = 
4T 3 = 4A 3 /g 2 , and the remaining half comes from the first term of the scalar potential 
in Eq.(|3|) containing derivative of the superpotential. The bulk cosmological constant 
Abuik = — 3T 6 /Mp comes from the last term — 3e^' 2 / M p|P| 2 /M| of the scalar potential 
in Eq.(|3|). These cosmological constants automatically satisfy the relation Aboundary = 
4M P a/— Abuik/3 which was implemented as a fine-tuning in the orbifold model 0. Our 
solutions (0) replace the boundary cosmological constant at a fixed point of the orbifold 
by a smooth physical wall configuration of the scalar field completely. The existence of the 
thin wall limit suggests that the M = 1 supergravity in four dimensions is consistent even 
on the orbifold, although we do not know if this point has been demonstrated similarly 
to the five- dimensional case ||. 

We can use the "effective supergravity" to extend our model to arbitrary spacetime 
dimensions D 0. Assuming the periodic potential W((j>) = (A D ~ 1 /g 2 ) sm(g(j)/ A^ D ~ 2 ^ 2 ) 
can be used, we obtain the D dimensional BPS wall solution connecting from = —tt/2 
to = tt/2 
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with Md as the Planck mass in D dimensions. A BPS wall different from ours has been 
obtained in a similar model in five-dimensional supergravity fTTfl. 



Classification and global properties of Non-BPS Solutions 



4 



In order to study possible multi-wall configurations, especially non-BPS configurations, 
we need to solve the equations of motion which are second order differential equations. 
It has been shown that these second order equations are equivalent to a set of first order 
differential equations which resemble the BPS equations provided there is only single 
scalar field forming walls |J . 

Let us assume that initial conditions for complex scalar 4>{y) , d<p(y) / dy and the warp 
factor A(y) , dA(y) / dy are real at some y, say y = — oo. The equations of motion for 
and A in our model then dictate that both and A are real for any y. Therefore we 
can ignore the imaginary part of the complex scalar field as long as we are interested in 
classical solutions. The bosonic part of the Lagrangian is given by 

We can now apply the method of Ref. || to our model taking as a real scalar field. Given 
the scalar potential (114]), we should find a real function W(<f>) by solving the following 
first order nonlinear differential equaiton 



( dW(<P) Y 
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V&)= — -772^(0)- (15) 



V d(f> J M, 



Then (f)(y) and A(y) are obtained by solving the following two first order differential 
equations 

dy dxf> • dy Ml W 1 ' 

It has been shown that these equations are equivalent to the equations of motion ||. 
It is convenient to rewrite these equations in terms of dimensionless variables 



dW 



±JV((f)) + 3a 2 W((j)) 2 , 1/(0) = cos 2 - 3a 2 sin 2 (18) 
d(j) 

d i = *m , d 4= -« 2 ^ ( 0). d9) 

dy d<p dy Zi 

To obtain a real solution, the right-hand side of Eq. (|18|) has to be real : V(<p) + 

3a 2 W 2 > 0. The shaded region in Fig.^ shows the forbidden region in W plane. It is 

easy to see that the exact (anti-)BPS solutions are those lines whose starting and ending 

points are both tangent at the top or bottom of these forbidden regions as shown in Fig.[| 

Any solution is represented by a curve in W plane. Since the first order differential 

equation ( |I8"D determines the solution uniquely once a value of W is specified at some 

as an initial condition, we observe: 1) No solution curves can intersect. Entire allowed 

region is covered by the curves once. 2) Curves can end only at the boundary of forbidden 

region or at infinity. 

Taking the positive sign for dW /dip in Eq.([H|), the asymptotic behavior of solution 
curves are given by W ~ ±exp(±v3o;0) — > ±oo. As illustrated in Fig.||, there are two 
different types of curves except the BPS solution curve : 
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Figure 2: Solution curves of Eq.(|l8|) with dW/d<p > in (f)W plane. Forbidden region is shaded. 

1. Starting at the boundary of forbidden region and ending at 0, W — > oo. 

2. Ending at the boundary of forbidden region and starting at 0, W — > — oo. 



Since the asymptotic behavior 0, W — > ±oo gives singularities || , |T2| at finite y (y ~ 
y ±ao =p e T< ^/(3a 2 ) as — > ±oo), we conclude that BPS wall solution ( |T2| ) is the only 
solution which are regular in the entire region of y if we do not allow any external sources 
for scalar as well as gravity, such as a boundary cosmological constant. If one looks 
into 0, W plane only, W = sin may appear a multi-wall winding solution. However, 
wall solution should be considered as a function of our base space — oo < y < oo. The 
solution curve approaches to the top or bottom of the forbidden region only at y — > ±oo. 
Therefore only a segment of W = sin between the top and the bottom of the forbidden 
region constitutes a wall solution which is a single wall. The non-existence of multi-wall 
BPS solutions is perhaps a peculiarity of our model which is in interesting contrast to the 
monopole or instanton case. 

We now wish to find a non-BPS solution by allowing a source due to a boundary 
cosmological constant XiS(y — yi) at appropriate points y = yi in the extra dimension into 
the Lagrangian fll4|) similarly to the orbifold model || J|]. We require that the equation 
of motion for the scalar field is satisfied without any external sources. Then only the 
equation of motion for the warp factor is modified to give the boundary condition 

dJ\. , . dA . . 1 < i 

^ + £ >-^- £ > = -2M? A <- < 20 > 



The boundary condition combined with Eq . (jig ) relates the boundary cosmological con- 
stant Aj at the connection point = 0» to a discontinuity of W at = 0j 

X l = 2(W(<p i + e)-W(<p i -e)). (21) 

Equation of motion for without external sources requires both and d<p/dy = dW /dip 
to be continuous across y. t . 

Let us denote the solution of Eq . (p~8|) specified by the initial condition (00,^0) as 
W±(0; (0o, Wo)) with the suffix +(— ) corresponding to the sign of dW/d(j). Symmetry of 
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Eq.(|lD gives 

B7 ± (0; (0o, Wo)) = W ± (0 + tt; (0 O + tt, W q )) = 



-W±(-0;(-0o,-W o )), (22) 



W+(0; (0o, Wo)) = W_(-0; (-00, Wo)) (23) 

As a simple solution to satisfy the requirement of continuity of and dip/dy = dW/d(p, we 
can choose to switch solution curves at = 7r from W(0; (0, W > 0)) to W(0; (2tt, —Wo)) 
and then at = 2tt to W(0; (27T, Wo)). As shown in Fig.|3] this solution is periodic in 
with the periodicity 2tt, and has a unit winding number and the symmetry S 1 /Z 2 which 
is consistent with the orbifold compactification as in Refs.|§, ||. We see that we cannot 
avoid negative cosmological constant (negative energy density). For small Wo, the solution 
consists of a BPS wall and an anti-BPS wall at fixed points = and = 7r of the orbifold 
supplemented by a small positive cosmological constant at = and a large negative 
cosmological constant at = tt. We can similarly construct solutions with arbitrary 
winding number with or without the orbifold symmetry by connecting solution curves at 
= utt with arbitrary integer n. These solutions may be regarded natural in the sense 
that fixed points occur where the scalar field is localized as walls. 




Figure 3: Non-BPS solution with unit winding number, (a) 4>W plane, (b) 0(y), and (c) A(y). 
Discontinuities in W(4>) at = 0, ±7r correspond to boundary cosmological constants. 




Figure 4: Non-BPS solution with unit windin; 
constant at = tt/2. (a) 0VP plane, and (b) 0( 
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If we allow the fixed points of orbifold at the mid point between two walls, we can 
connect solution curves at = mr + tt/2. For instance we can switch solution curves at 
= tt/2 from W(<j>\ (0, W > 0)) to W$\ (vr, -W )) and then at = vr to W(<j>\ (vr, W )). 
This configuration has a periodicity = + tt and has a symmetry of S 1 /(Z 2 X Z 2 ) 
orbifold as shown in FigfJ. In the limit of small H 7 ^ this solution can be regarded as a 
model essentially replacing the positive cosmological constant by a smooth physical wall 
configuration of scalar field whereas the negative cosmological constant remains almost 
the same as in the original model 0,0]- Actually we can also connect solution curves at 
arbitrary points in by switching to the curve with the same derivative at the same 0. 

Another possibility to connect different solution curves is to switch to curves with 
opposite sign dW/dtfi. This is possible only at the boundary of the forbidden region 
where dW /dip = 0. If we connect solutions at the boundary of the forbidden region, we 
can do without putting boundary cosmological constant by switching to the curve starting 
from the same value of W with the opposite sign for dW/d<ft. A typical example with no 
winding is shown in Fig.[|(a). This solution also has the orbifold symmetry S l /Z 2 . 




Figure 5: Non-BPS solutions without winding number in <j>W plane. Discontinuities in W(cf>) 
at = 0, ±7r correspond to boundary cosmological constants, (a) Smooth with two types 
of boundary cosmological constants at y = 0, ttR/2. (b) with only one type of boundary 
cosmological constant but with a break of scalar field at y = nR/2. (c) Using a piece of the BPS 
solution and boundary source for 0. 



If we allow boundary cosmological constants at the boundary of the forbidden region, 
we can connect solutions by switching to the curve starting from the W as well as dW /dcf) 
of opposite sign at the boundary of forbidden region. We illustrate the simplest of such 
solution in Fig.|5|(b). This solution has the symmetry S 1 /(Z 2 x Z 2 ) and has no winding 
number. If we relax our principle and allow the boundary action to contain a source term 
for the scalar field 0, we have more varieties of solutions. For instance we can even use 
the BPS solution to construct a solution, simplest of which is illustrated in Fig.|5](c). It 
has the orbifold symmetry S 1 /Z 2 but without any cosmological constant at one of the 
fixed point where the scalar field energy density is localized. This is a model replacing 
the positive cosmological constant of the original model || completely by a smooth thick 
wall configuration of scalar field at the cost of having a boundary source term for the 
scalar field in the other fixed point. We can also construct more complicated solutions 
with or without the orbifold symmetry by combining these various types of boundary 
cosmological constants. 
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The precise amount of the cosmological constants and the distance between the walls 
are determined by Eqs. (|18"D and fl20P . We will give a systematic analytic approximation 
valid for small gravitational coupling below. 

A systematic approximation for weak gravity 

To establish the limit of no gravity (a — > 0) for the first order differential equations 
( |I8D and (|19D, let us first examine the equation of motion for scalar field without gravity 
using the scalar potential V(<p; a = 0) in the absence of gravity 



dV(4>; a = 0) _ d?ip 
dy 



2-f 2 . (24) 



By multiplying d<p/dy, we can integrate it once with an integration constant wq 

V(4>;a = 0)= -3w 2 . (25) 



This first order differential equation can be identified as the no gravity limit of Eq. fllq) , 
if a function W(<f>) is defined by dW/dcj) = d<p/dy. Because of Eq . (|i"8|) , Eq . (|T9|) dictates 
that the limit of vanishing gravitational coupling is given by 

a -»• 0, aW(4> = 0)=w fixed. (26) 

Eqs.(^) and flT8j) show that wo = corresponds to the BPS solution. 

Since Wo = W((j) — 0) ~ 1/a — > oo, we define a function W more appropriate for an 
expansion in powers of gravitational coupling a 



n=0 n=0 n=0 

The first order differential equations appropriate for an expansion in powers of a is 



2^ = ^ = ±Jv(4>) + 3(w + aW(4>)) 2 , ^ = -aw - a 2 W(4>). (28) 
d(f) dy V V J dy 

It is interesting to note that the first correction is of order a = A/ (gMp) rather than the 
usual a 2 . For our specific model of gravitationally corrected periodic potential, we obtain 

^ = ^ = J^Ti< ^ = , - (29) 

d( t> ^ V # d » y/cos^ + 3^ 2 

-r^ = 0, — - = -au> , — = -qV°'W,---. 30 
dy dy dy 

The zero-th order equation for (p^°\y) in ( p9|) gives precisely our exact non-BPS solution 
which was obtained in the case of global SUSY [17j 



^\y) = ^rn y —^,k), k= J_ , (31) 
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W (0 \4>) = \ f d9Vl-k 2 sin 2 9 = jE U, k) , (32) 
k Jo k \ / 

where am is the amplitude function and E is the elliptic integral of the second kind. 
The warp factor to the zero-th order is given by an irrelevant normalization A^\y) = 
constant. 

Using Eq. (|21|) , we can now determine to the lowest order in weak gravitational coupling 



the boundary cosmological constant needed to satisfy the boundary condition (|20|) . For 
instance, the unit winding solution in Fig.^| requires the boundary cosmological constants 
Ao at cf) = and Ai at cf) = it 

M P A 2 , M P A 2 A 3 . 

A = 4^— w , Ai = —4 — - — wo - S—E(k), 33 
9 9 9 2 k 

where E(k) is the complete elliptic integral of the second kind E(k) = E(tt/2, k). Similarly 
other situations can also be expressed in terms of elliptic integrals. 

A systematic approximation in powers of gravitational coupling is given by Eq.(|29[). 
For instance the first order correction is given by integrating the second of Eq. (p9|) 

fi - ft = r d) ilw„mk)' (34) 

W® = (* M IpEML, A«\y) = - ama y. (35) 
Jo VI - k 2 sin 2 6 

We can obtain gravitational corrections to any desired order with increasing complexity. 
We stress that SUSY is not needed for our method to work similarly to Ref . || . Therefore 
it should be useful to obtain solutions for theories with or without gravity irrespective of 
presence or absence of SUSY. 
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Figure 6: Solution curves of Eq.(|l8[) in (f>W plane with larger gravitational coupling a 
A/(gMp) = 0.01, • • • , 1.0. The curve stays almost the same for a < 0.01. 



Let us also examine the situation with stronger gravitational coupling. The curve 
stays almost the same for a < 0.01. We see in Fig.^ that the larger values a > 0.4 
reveal an early emergence of singularities even for relatively smaller values of <fi. However, 
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qualitative features remain quite similar for small values of 0. Therefore we expect non- 
BPS multi-wall configurations are still possible for medium gravitational coupling. It is 
interesting to note that our BPS solution exists even with strong gravitational couplings 
at least classically. 

Finally we wish to note that the existence of the global SUSY limit fl3lD , fl32|) implies 
that we are sure that there is a mass gap for fluctuation modes of the scalar fields. There 
is a possibility for our wall scalar field <fi to serve as a Goldberger-Wise type stabilizer 
[PH. We will settle the issue of stability by analyzing the possible new modes associated 
with the gravity supermultiplet in a subsequent report. 
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